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We develop models for the propagation of intense pulses in solid state media which can have 
either saturated absorption or exhibit reverse absorption . We show that the experiments of Bigelow 
et a/.[Phys. Rev. Lett. 90, 113903 (2003); Science 301, 200 (2003).] on subluminal propagation in 
Ruby, and superluminal propagation in Alexandrite are well explained by modelling them as three 
level and four level systems respectively, coupled to Maxwell equations. We present results well 
beyond the traditional pump-probe approach. 



Since the discovery of ultraslow light with a group velocity 17 m/sec in a Bose condensate by Hau et al. [1-3] many 
experiments have reported slow light in a variety of media [4-6] . Kash et al. [4] demonstrated light propagation with 
a group velocity of 90 m/sec at room temperature in Rb vapor. Using Zeeman coherences, Budker et al. [5] reported 
slow light with group velocity 8 m/sec in Rb vapor. Hemmer et al. [6] first reported slow light in solid state material 
namely Pr doped Y2SiC>5, maintained at a cryogenic temperature of 5K. The slow light ideas have been successfully 
used in storage and retrieval of light pulses [7,8]. The understanding of storage and retrieval of light pulses has been 
provided by Dey and Agarwal [9], using the adiabaton theory of Grobe, Hioe and Eberly [10]. 

Work on pulse propagation continues to produce interesting results [11-14]. Recently, Bigelow et al. [15] showed 
the propagation of light pulses in Ruby at a group velocity of 57.5 m/sec. This experiment differs considerably from 
all earlier experiments which were based on clcctromagnetically induced transparency [16]. Bigelow et al. recognize 
that a two level system driven by a strong field and a probe gives rise to a hole in the probe response function with 
a width of the order of 1/Ti, where T\ is the longitudinal relaxation time [17]. Note that in a material like Ruby 
the transverse relaxation time X2 <C T%. These authors also discovered that they need not use separate pump and 
probe fields. A field with peak power of the order of saturation intensity could be slowed down considerably to about 
57.5 m/sec. Further Bigelow et al. [18] extended their work to a material like Alexandrite, where they reported 
superluminal propagation. There have been several earlier reports of superluminal propagation in solid state systems 
[19]; and in vapors [20,21]. 

The purpose of this letter is to study the propagation of intense pulses in a homogeneously broadened medium, like 
Ruby, which can exhibit saturated absorption or a medium like Alexandrite, which can exhibit reverse absorption. 
Note that the pulse propagation in a nonlinear transparent medium has been extensively studied [22]. The systems 
studied here differ from the previous studies as our systems posses very strong transverse and longitudinal relaxation 
effects. In order to model the experiments, we model Ruby as a three level system and Alexandrite as a four level 
system. We solve the system of coupled equations numerically to delineate various aspects of pulse propagation. We 
do not make any approximation on the strength of the pulses so that we can model experimental observations on 
strong pulses. We calculate group velocity from the relative delay or advancement between the reference pulse and 
the output pulse. We present numerical results on the propagation of Gaussian and modulated pulses and show good 
agreement with the experimental data of Refs. [15,18]. The experiments of Bigelow et al. fall in two categories. One 
consists of a weak probe and a strong coherent cw field. These were explained in terms of the response to a probe 
field of a two level medium pumped by a coherent field [23] . The other category consists of the self delay of a strong 
pulse. The latter requires solutions of the coupled Maxwell-Bloch equations and this is the case we concentrate on. 

Media with saturated absorption — -For pulse propagation in the three level model in Ruby, we represent the ground 
state as \g), the 4 i 7 2 absorption band as |ei) and the levels 2A and E as \e2)- In Ruby one has very rapid decay of the 
level |ei) to |e2), and as a result some of the coherences become irrelevant on experimental time scale. The density 
matrix equations for the model of the Fig. 1 are 
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where p., = (e ( |p|e 3 .), i,j = 1,2. The Rabi frequency 20 is defined by 2Q(z,t) = 2d lg ■ £(z,t)/h, where d lg is the 
dipole matrix element and £{z, t) is the envelop of the pulse. We assume that the carrier frequency, to, is on resonance 
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with the frequency of the |ei) < — ► |g) transition. Under the approximations, Ti 3> r 2 ,f2; p lg ~ 0, we derive the 
approximate equation for the evolution of the ground state population as 

p ga = 2Y 2 (l- Pgg ) - 2^p gg (2) 
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Note that we can prove that p lx w 0, if Ti r 2 ,0. Under the same conditions and the slowly varying envelop 
approximation, the evolution equation for the Rabi frequency of the field is governed by 



dCl a ^ 
~dz~~~^2 



Clpgg, Cl = Cl/Cl sat , (3) 



where an — 4iru\d lg \ 2 /cHTi and Cl sat = 2^TiT 2 . In Eqs. (2) and (3) we have used the pulse coordinates i.e, t — z/c, z. 
The time derivative in Eq. (2) is with respect to (t — z/c). The time t can be expressed in units of 1/2IV For 
numerical computation, we consider two types of input pulses, viz, a Gaussian pulse with a temporal width > l/r 2 

ci in = n° e [-* 2 / 2 - 2 ] (4) 

and amplitude modulated pulse 

n?„(t)=I = I (l + mcos[At]). (5) 

The Equations (2) - (5) are our working equations. We use these for numerical computations. We calculate the 
evolution of the pulse for arbitrary values of Cl or Io- Some typical results for the Gaussian pulses are shown in the 
Fig. 2. We get group velocities in the range 50 m/sec for Cl/Cl sat ^ 1 and the transmission is rather small. In Fig. 
3, we exhibit the behavior of v g and transmission as a function of the input intensity. These results are in agreement 
with the experimental findings of transmission in the range 0.1%. In Fig. 3, we also show for comparison the results of 
the group velocity and the transmission for the propagation of an intense pulse through a two level system described 
by the traditional Bloch equations. The coupled Maxwell-Bloch equations under the approximation T\ 3> T 2 arc 
given by 

PiiT 1= - Pll + 2\Cl\ 2 (l-2 Pll ) (6) 
P 99 + P 11 = 1 (?) 
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where Cl = Cl^/T\T 2 and dot denotes d/d(t — z/c). As seen from the Fig. 3, there are substantial differences in the 
propagation of pulses in two level and three level media. Note that the time T\ is equal to l/2r 2 . We believe that, 
in the light of the energy level diagram of Ruby, it is more appropriate to model it as a three level system. 

We next consider input pulse as a modulated pulse given in Eq. (5). The output pulse is modulated with a phase 
shift (time delay). We show this time delay as a function of modulation frequency for two different pump powers in 
the Fig. 4. The results in Figs. 2-4 are in excellent agreement with the experimental data(c/for example Fig. 3 of 
Rcf. [15] with our Fig. 4). 

Media with reverse saturation — For the supcrluminal propagation in Alexandrite, Bigelow et al. recognized how 
the reverse saturation mechanism [24] can be at work in a material like BeAl 2 04 doped with Cr 3+ ions and with 
some Cr 3+ ions replaced by Al 3+ . The reverse saturation produces an antihole in the susceptibility for the probe in 
presence of a pump field. The antihole can result in the supcrluminal propagation. In what follows we show how the 
measurement can follow by modelling the system as a four level system to account for reverse absorption. The model 
is shown in the Fig. 5, where state 4 A 2 as \g), the absorption bands 4 T 2 and 4 Ti as |ei) and the level 2 E as |e 2 ). The 
density matrix equations are now given by 
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Following the same procedure as in the case of Ruby, we have derived the working equations 
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where ao gives the reverse saturation. Following the experimental data of Bigelow et al. [18], wc estimate (ao/ao) ~ 4 
. The Eqs. (10) and (11) are numerically integrated for the input Gaussian pulse given by Eq. (4). A representative 
set of results is shown in the Fig. 6. This Figure also shows how the group velocity and net transmission depends 
on the peak intensity of the Gaussian pulses. It should be borne in mind that in the range of the intensities of Fig. 
6, no perturbation theory can be used. One has to study the full nonlinear behavior. We also notice that the input 
pulses get distorted in shape. The distortion becomes more pronounced as the nonlinearity of the medium becomes 
more pronounced. 

In conclusion, we have shown how to model the propagation of intense pulse in solid state media with very strong 
relaxation effects. The media can exhibit either saturated absorption or reverse absorption. Our modelling goes well 
beyond the traditional pump-probe approach. We specifically present results on the propagation of pulses in Ruby 
and Alexandrite. Our model would also be applicable to other systems where reverse absorption could be dominant. 
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FIG. 1. Three level model for Ruby Crystal 

FIG. 2. The solid curve shows light pulse propagating at speed c through 7.25 cm in vacuum. The long dashed and dot-dashed 
curves show light pulses propagating through a medium of length 7.25 cm at different input amplitudes. The temporal width 
a of the Gaussian pulse is 20 msec and 1/2F2 = 4.45 msec. The part (b) gives the amplitudes, of the output pulse normalized 
to the input amplitudes. The transmission increases with increasing the input field intensity. 

FIG. 3. Variation of transmissions and group Velocities as functions of the input amplitude of the light pulse. The solid 
(dashed) curve gives the intensity transmission of the pulse for the two (three) level model of the medium. The corresponding 
group velocities are given by the dotted curve (two level model) and the long dashed curve (three level model) . The light pulse 
is propagating through the medium of length L=7.25 cm. 

FIG. 4. Time delay of the light pulse as a function of modulation frequency for two different input powers. The modulation 
index, m, is equal to 0.05. Note that the output pulse can be fitted to the from I(t) = I (t)(l + mcos(At + 0)). The overall 
transmission is small, e.g., I (t)/I s = 0.00165 for A/2F2 = 2 



FIG. 5. Four level model for Alexandrite crystal 

FIG. 6. The solid curve of (a) shows light pulse propagating at speed c through a distance of 7.25 cm in vacuum. The 
dotted, long dashed and dot-dashed curves depict light pulse propagating through a medium of length 7.25 cm at different 
input amplitudes. The pulse width a is 500 pisec, whereas l/2r2=250/isec. Fig (b) shows the amplitude of the output pulse 
normalized with input amplitude. The transmission is decreased on increasing the input field intensity . 
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